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D ■ Abstract 



We construct a well-defined lattice-regularized quantum theory formulated in terms of funda- 
mental fermion and gauge fields, the same type of degrees of freedom as in the Standard Model. 
The theory is explicitly invariant under local Lorentz transformations and, in the continuum limit, 



P-i under diffeomorphisms. It is suitable for describing large nonperturbative and fast-varying fluc- 



tuations of metrics. Although the quantum curved space turns out to be on the average flat and 



CN " smooth owing to the non-compressibility of the fundamental fermions, the low-energy Einstein 

>[ 

limit is not automatic: one needs to ensure that composite metrics fluctuations propagate to long 



distances as compared to the lattice spacing. One way to guarantee this is to stay at a phase 



OO " transition. 

o; 

We develop a lattice mean field method and find that the theory typically has several phases 
in the space of the dimensionless coupling constants, separated by the 2 nd order phase transition 
surface. For example, there is a phase with a spontaneous breaking of chiral symmetry. The effec- 
tive low-energy Lagrangian for the ensuing Goldstone field is explicitly diffeomorphism- invariant. 
We expect that the Einstein gravitation is achieved at the phase transition. A bonus is that the 
cosmological constant is probably automatically zero. 
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I. INTRODUCTION 



We live in a world with fermions, and they must be included into General Relativity. 

fir 

The standard way one couples Dirac fermions to gravity is via the Fock-Weyl action |2J: 
Fermions interact with the frame field (also known as vierbein, repere or tetrad) and with 
the spin connection co^ B being the gauge field of the local Lorentz group. The frame and 
the spin connection are a priori independent field variables. The bosonic part of the action 
has to be written through e M and accordingly. This is known for the last 90 years as 
Cartan's formulation of General Relativity {3]. Speaking generally, it is distinct from the 
classic Einstein-Hilbert formulation based on the Riemann geometry, since it allows for a 
nonzero torsion. We stress that the presence of fermions in Nature forces us to make a 
definite choice in favor of the Cartan, as contrasted to the Riemann geometry. 

In practice, however, it is hardly possible to detect the difference. In the leading order in 
the gradient expansion of the gravitational action written down in terms of e M and u^, the 
saddle-point equation for says that torsion is on the average zero. Therefore, Cartan's 
theory reduces to that of Einstein. 

In the next order in p 2 /Mp where Mp is the Planck mass and p is the characteristic 
momentum, a four-fermion contact interaction appears from integrating out torsion. Its 
strength is many orders of magnitude less than that of weak interactions J4] therefore this 
correction will hardly be detected any time soon in the laboratory. In principle, it modifies 
e.g. the Friedman cosmological evolution equation that follows from the purely Riemannian 
approach. However the correction remains tiny as long as fermions in the Universe have 
Fermi momentum or temperature that are much less than Mp If they reach that scale 
such that the four-fermion correction becomes of the order of the leading stress-energy term, 
the theory itself fails since the gradient expansion (5), (| from where it has been derived, 
becomes inapplicable. There is no agreed upon idea how the theory looks like at the Planck 
scale; in particular, quantum gravity effects are supposed to set up there. 

Being indistinguishable from Einstein's equation in the range where observations are 
performed, Cartan's theory, however, has a critical feature when one attempts to quantize 
it. The bosonic part of the action is written in terms of e M and u^. To preserve the 
required general covariance or invariance under the change of coordinate system, called 
diffeomorphism, any action term is necessarily odd in the antisymmetric Levi-Civita tensor 
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e K\nv ^hat m akes all possible diffeomorphism-invariant action terms not sign-definite [7J. 

The simplest example is the invariant volume itself or the cosmological term, J d A x det(e). 
If the frame field is allowed to fluctuate, as supposed in quantum gravity, the sign of det(e) 
can continuously change from positive to negative or vice versa. Of course, det(e) = is a 
singularity where the curved space effectively looses one dimension but it is not possible to 
forbid such local happenings in the world with a fluctuating metric, see the illustration in 
Fig. 1, left. Moreover, if det(e) goes to zero linearly in some parameter t, it has to change 
sign by continuity, see Fig. 1, right. The same is true for any diffeomorphism-invariant 
action term. 




det(e)>0 det(e)<0 det(e) > 1 e \ 

det(e) > det(e) = det(e) < 

FIG. 1: Left: An example of a space with alternating sign of det(e); Right: det(e) changes sign by 
continuity of the frame field. 

In the standard Riemannian formulation, one writes the invariant integration measure 
with the help of ^det(g) where g^ u = e^e^ is the metric tensor, hence det(g) = (det(e)) 2 is 
sign-definite. Its square root, however, should be understood as y / detTo) = det(e) and can 
have any sign. If it passes through zero it changes sign by continuity jsj. 

This fundamental pathology of any diffeomorphism-invariant quantum theory has not 
been stressed before, probably for two reasons. First, one commonly deals with the pertur- 
bative quantization about flat or e.g. de Sitter metric such that the main concern is the 
absence of runaway fluctuations from that point only. However, when quantizing gravity, 
one has to be concerned with large non-perturbative fluctuations as well. Second, usually 
Minkowski space-times are considered where the integration measure exp(i Action) is oscil- 
lating anyway independently of the action sign. However, a theory with a sign-indefinite 
action in Euclidian space where the weight is exp(— Action) is usually fundamentally sick 
also in Minkowski space. An illustration is provided by the scalar <p 3 theory, see Fig. 2. 
Perturbation theory exists there in the usual sense near = 0. However, if in Euclid- 
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ian space the theory does not exist, in Minkowski space one cannot define properly the 
non-perturbative Feynman propagator. There will be also other pathologies related to the 
possibility of tunneling to a bottomless state. 




FIG. 2: The </> 3 theory is fundamentally sick both in Euclidean space where it is unbounded, and 
in Minkowski space where it can tunnel to a bottomless state. 

In gravity theory, Euclidian formulation has its own right, for example in problems related 
to thermodynamics and to tunneling, like in the Hawking radiation problem where paradoxes 
are encountered just because we do not know how to quantize Euclidian gravity. If a theory is 
well defined for Euclidian signature, it is usually possible to Wick-rotate it to the Minkowski 
world. Therefore, for clearness we shall discuss here Euclidian gravity. 

Any diffeomorphism-invariant action, with any number of derivatives, is not sign-definite 
in Euclidian space and hence cannot serve to define quantum gravity non-perturbatively. 

At this time, we see only one way to overcome the sign problem, and that is to use in part 
fermionic variables in formulating quantum gravity microscopically, rather than only bosonic 
ones. Integrals over anticommuting Grassmann variables are well defined irrespectively of 
the overall sign in the exponent of a fermionic action. The reason is that in fermionic 
integrals introduced by Berezin 9j one actually picks up only certain finite order in the 
Taylor expansion of the exponent of the action, such that the overall sign does not matter. 



One calls it spinor quantum gravity. It has 
recently in a series of papers by Wetterich 



3een advocated by Akama lCj , Volovik [ll|] and 
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15 



on other grounds. 

lo) ]) that at the fundamental, microscopic 



More specifically, we suggest [?J (see also Ref. 
level gravity theory is defined as a theory of certain fundamental anticommuting spinor fields 
We wish to preserve local gauge Lorentz symmetry exactly at all stages, and for that 
we need the explicit connection field u^. The frame field and the metric tensor g^ will 
be composite fields making sense only at low energies. The basic independent variables will 
be ip\ip and the gauge field u^, the same type of degrees of freedom as in the Standard 



4 



Model. We believe that using the same type of variables as in the Standard Model will help 
to unify all interactions [7|. As far as only gravity is concerned, the fundamental spinor 
fields m ay or may not be related to the fundamental matter fields. We introduce the 
main building blocks of the theory in Section II. 

A quantum field theory is well defined if it is regularized in the ultraviolet. We shall 
regularize spinor quantum gravity by introducing simplicial lattice (made of triangles in 
2 dimensions, tetrahedra in 3 dimensions, 5-cells or pentachorons in 4 dimensions, etc.) 
covering an abstract space, such that the simplex vertices are characterized and counted by 
integers i. Only the topology of this abstract number space matters, e.g. the number of 
nearest neighbors, etc. 

Each vertex i in the number space corresponds to the real world coordinate by a cer- 
tain map x^(i). Diffeomorphism-invariance means that the theory should not depend on 
the coordinates x^(i) we ascribe to the vertices. Also, in the continuum limit (implying 
slowly varying fields) the action should be of the form J d d x C(x) and invariant under dif- 
feomorphisms, x M — > x l)l (x). The integration measure over the fields in the path integral 
formulation should be also diffeomorphism-invariant. In addition, we require exact gauge 
invariance under local Lorentz transformations. We build fermionic actions satisfying these 
conditions in Section III, and regularize them by putting on a lattice in Section IV. 

After constructing a completely well-defined lattice-regularized quantum theory, the next 
question to address is whether the continuum limit can be achieved and whether it reduces to 
the Einstein-Cartan theory in the low-energy limit. The continuum limit is obtained when 
and if field correlations spread over a large distance in lattice units. The trouble is that the 
quantum theory one deals with in this approach is a typical strong-coupling theory where 
most of the correlations die out over a few lattice cells. Contrary to the standard lattice 
gauge theory where long-range correlations are ensured by simply taking the weak-coupling 
limit /3 — > oo, in spinor quantum gravity there is no obvious handle to make the correlations 
long-ranged. 

The main trick and the invention of this paper is to ensure long-range correlations by 
adjusting the bare dimensionless coupling constants to the point (or line, or surface) where 
the theory undergoes a phase transition of the 2 nd kind. At such point, all correlation func- 
tions become long-ranged, and the Einstein theory will be guaranteed in the low-momenta 
limit by the inherent diffeomorphism invariance. 
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Second-order phase transitions occur in theories where there is an order parameter, usu- 
ally related to the spontaneous breaking of a continuous symmetry. Our primary goal in 
this project is to demonstrate that second-order phase transitions are typical in the kind of 
diffeomorphism-invariant theories we consider. We develop in Section V an original mean 
field method well suited for the search of the phase transitions, and check its accuracy in 
the Appendix where it is probed in an exactly solvable model, with very satisfactory results. 
Using this method we unveil the phase diagram of a generic 2-dimensional lattice spinor 
gravity, in the space of the bare coupling constants, Section VI. The model has two contin- 
uous symmetries: the U(l) chiral symmetry and the U(l) symmetry related to the fermion 
number conservation; both can be in principle spontaneously broken. 

It turns out that there is a range of bare couplings where the fermionic lattice system 
experiences spontaneous breaking of chiral symmetry. In the particular model we studied we 
did not observe spontaneous breaking of the fermion number, however it can happen in other 
models. This is an interesting finding per se but it may be also of use in the attempts to 
unify quantum gravity with the Standard Model. On this route, one expects one or several 
spontaneous breakups of continuous symmetries. 

The 2-dimensional model we consider in some detail has certain nice features. First, 
the physical (invariant) volume (V) is extensive i.e. proportional to the number of lattice 
points taken. This is not altogether trivial since nonperturbative metric fluctuations allow, 
in principle, "crumpling" of the space, and that is what some researchers indeed typically 
observe in alternative nonperturbative approaches to gravity. In spinor gravity, it is a natural 
result following from the non-compressibility of fermions. Second, the quantum average of 
the curvature turns out to be zero such that the empty space without sources is effectively 
flat. This is also a welcome feature since the natural result in nonperturbative gravity is that 
the curvature is of the order of the cutoff, that is of the Planck mass, which is unacceptable. 
Third, despite flatness the theory definitely describes a fluctuating quantum vacuum, as 
exemplified by the fact that the physical volume variance or susceptibility (V 2 ) — (V) 2 is 
nonzero. 

As a result of the spontaneous breaking of continuous symmetry (here: chiral symmetry), 
a Goldstone field appears. We check by an explicit calculation in Section VII that the 
low-momentum effective ("chiral") Lagrangian for the Goldstone field is diffeomorphism- 
invariant as expected. This invariance is rooted in the way we construct the original lattice 
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action for spinors. The appearance of a Goldstone particle means that a definite bilinear 
combination of fermions is capable of propagating to large distances. However this is not 
enough: in order for the system to totally loose memory about the original lattice, all degrees 
of freedom have to propagate to long distances in lattice units. This happens only at a phase 
transition where we expect that the Einstein-Hilbert action emerges as a low-energy effective 
action for the classical metric, with the cosmological constant being automatically zero, see 
Section VIII. In Section IX we discuss the dimensions of various quantities and fields used 
throughout the paper. We summarize in Section X. 



II. COMPOSITE FRAME FIELDS 

Following Ref. [3] we introduce a composite frame field built as a bilinear fermion 
"current". In d dimensions the frame field transforms as a vector of the SO(d) Lorentz 
gauge group: 

e» 0**(x)e*(x). (1) 

Since A,B, . . . = 1, . . . , d are flat Euclidean indices we can equivalently write them either as 
subscripts or superscripts. The frame field transforms also as rank-one tensor (world vector) 
with respect to diffeomorphisms x M — > x' M (x): 



efa) diffcomorphism ) (2) 



Our basic objects are fermion fields ip(x), x/j'(x) assumed to be world scalars under diffeo- 
morphisms, and transforming according to the spinor representation of the Lorentz group, 

ip(x) ->• V{x)ij}{x), ifl{x) ->• tf{x)V\x), V G SO(d), 

ip, ip\x) — > ip, ijj^ (x'(x)) . (3) 

n 

The dimension of the spinor representation is df — 2' d//2 l, see, e. q.[17j]. 
We introduce the covariant derivative in the spinor representation, 

V, = 8, - l ~ufX AB , = \< B ^AB (4) 

where u AB is the spin connection in the adjoint representation of the SO(d) group, and E^b 
are its df x df generators, 

Eab = ^[iaIb], (5) 



built from Dirac matrices 7^ satisfying the Clifford algebra, 

{1a1b] = 25 A b Idfxdf (6) 

In the adjoint (antisymmetric tensor) representation the corresponding covariant derivative 
is 

Df = d,8 AB + uf. (7) 
Its commutator defines the curvature 

[D„ D V \ AB = T AB , [V„ V,] = -ffiVAB, (8) 

where 

ft? = + ufdS* -(fx^u). (9) 

One can built two distinct bilinear combinations of the fermion fields, transforming as 
the frame field (jT121> : 

e A = 1(^ + ^7^), (10) 
f A = V>VV^ - ^T^- (11) 

To check that e A and f A transform as a vector (OQ) one needs the relation between the matrix 
V rotating spinors ([3]) and the matrix O rotating vectors, 

q ab = 1 Tr (ytyVy*) . 
d f 

Given that anticommute, the above bilinear operators are Hermitian. 

We can define the bilinear fermion operator that plays the role of the torsion field, for 
example, 

7>) ^ \ (D AB e B - D AB e B ) = \t ab (12) 
and similarly for the other composite frame field / M fllip . 

III. DIFFEOMORPHISM-INVARIANT ACTION TERMS 

One can now construct a sequence of many-fermion actions that are invariant under local 
Lorentz transformations and also diffeomorphism-invariant, using either e A or f A (or both) 
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as building blocks: 

k = 0,1,..., [d/2], 

where e WM2 '" w is the totally antisymmetric (Levi-Civita) tensor. Notice that So is the analog 
of the cosmological term but there are many of them since one can replace any number of 
e^'s by f A, s, Si is the analog of the Einstein-Hilbert-Cartan action linear in curvature, and 
the last action term SW21 for even d is a full derivative. Apart from full derivatives, there 
are 3 possible action terms in 2d, 6 terms in 3d, 8 terms in Ad, 12 terms in 5d, etc. 

The use of e MlM2 ' Atd is obligatory to support diffeomorphism-invariance. In principle, one 
can construct Lorentz-invariant action terms by contracting the flat indices with Kronecker 
deltas instead of e x " d , however that will make the action term P- and T-odd. For example, 
there is a well known P, T-odd term in four dimensions, called sometimes the Hoist action, 



kX/ii/ -pAB „A„B 



19 



20], but we do not consider such terms here. 
One can add to the list of admissible action terms any of the actions fflBl multiplied by 
any power of the world and Lorentz-group scalar (ip^ip)', we shall consider such kind of terms 
later on in relation to the spontaneous breaking of chiral symmetry. 

All action terms (1131) are apparently invariant under two global U(l) rotations: 

• phase rotation related to the fermion number conservation, ip — y e^ip, ^ — y i$e~ %OL 

• chiral rotation for even dimensions d, ip — y e ll3,yd+1 ip, ip^ ~^ w e ll3ld+1 , where 7^+1 = 
i d/2 lil2 ...id, {ld+ild} = 0, 7j +1 = 1, 

since both e A and f A are invariant under these transformations. The two corresponding 
Nother currents are conserved. However, both symmetries can be spontaneously broken by 
interactions, and we shall see that this is what indeed typically happens. 

IV. SPINOR GRAVITY ON THE LATTICE 

In order to formulate quantum theory properly one has to regularize it at short distances. 
The most clear-cut regularization is by (lattice) discretization, however diffeomorphism- 



invariance imposes severe restrictions on it, see recent discussion by Wetterich 13Ml5l|. We 
impose two basic requirements: 
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• explicit invariance under local gauge transformations of the Lorentz group, small or 
large (as in lattice gauge theory) 

• if the fields vary slowly in lattice units, i.e. in the continuum limit, the lattice action 
reduces to one of the diffeomorphism-invariant action terms ffTS]) and the like. 

A. Triangulation by simplices 

To that end, we introduce an abstract discretized space where only the topology of ver- 
tices and edges connecting neighbor vertices is chosen beforehand and fixed. We find that 
the simplest hypercubic topology does not work. Only in two dimensions it is possible, for 
accidental reasons, to fulfill item 2 above by introducing a square lattice. In higher dimen- 
sions, the simplest but sufficient construction is to use a simplicial lattice. For uniformity, in 
two dimensions we also consider a triangle lattice made of three- vertex cells. In 3d simplices 
are tetrahedra or 4-cells, in 4d these are pentachorons or 5-cells, and so on. 

It is always possible to cover the whole d- dimensional space by (d+ l)-cells or simplices, 
although the number of edges entering one vertex may not be the same for all vertices. 
Alternatively, the number of edges coming from all vertices is the same but then the edges 
lengths may vary, if one attempts to force the lattice into flat space. Since only the topology 
of the nearest neighbors matters and the abstract "number" space does not need to be flat, 
this is also acceptable. The important thing is that the chosen set of cells should fill in the 
space without holes and without overlapping [l8|. 

All vertices in a simplicial lattice can be characterized by a set of d integers. For brevity 
we label these d numbers by a single integer i. Each vertex has its unique integer label i, 
supplemented with a rule what labels are ascribed to the neighbor vertices forming elemen- 
tary cells. We shall denote the d + 1 labels belonging to one cell by i — 0, 1, . . . , d. In this 
section, we write down the full lattice action as a sum over actions for individual simplicial 
cells, therefore we shall not be concerned with the precise geometric arrangement of the 
cells. 

Each vertex in the abstract number space corresponds to the real world coordinate by a 
certain map x^{€). The goal is to write possible action terms in such a way that if the fields 
vary slowly from one vertex (or link) to the topologically neighbor one, the action reduces 
to one of the possible diffeomorphism- and local Lorentz- invariant action term in Eq. ffl3|) . 
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We start by writing the volume of an elementary cell (simplex) in a given coordinate 
system in d dimensions. It can be presented as a determinant of a d x d matrix, 

Simplex = ^det (Atii )(xf - a#), (14) 

where Xq is the coordinate ascribed to one of the vertices, and xf , i = 1 . . . d are the 
coordinates ascribed to all the other vertices. We introduce the notion of a "positive order" 
of vertices % in the cell: it is such that for smooth functions x% the volume (|14p is positive. 
An odd permutation of vertices in this set makes a "negative order" . 
It will be convenient to use the antisymmetric symbol 

if if. does not belong to a given cell, 

1 if the set i ,ii...i d is in the positive order, (15) 
— 1 if the set i ,ii...i d is in the negative order. 

With the help of this symbol the cell volume (fl4|) can be written as 

Simplex = ^ -q, — x i ) \ X i 2 — x i ) ■ ■ ■ ( x ij — x io) ■ (16) 

B. Lattice action 

The building blocks of our construction are anticommuting spinor fields that are 

world scalars and "live" on lattice vertices i, and the parallel transporter Uij. As in any 
lattice gauge theory, we replace the connection by a unitary matrix "living" on lattice 
links |3], 

Uij = Pexp £ ufH AB dx^ , V* = U} r (17) 

In terms of these lattice variables the discretized versions of the composite frame fields 
(flUyni) are: 



e& = i^)U 3 a A U l3 ^ - ^T^i), (18) 



ft = ^VMj - ^Pa^i- (19) 

The difference between e and / is that the first has both fermions in the same vertex whereas 
in the second fermions are residing in the neighbor vertices. 
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Expanding all fields in Eqs. (j!8|19p around the center of a cell x = ^-j- Yui=o x * we obtain: 

~4, = K-xf)e» + 0(Ax 2 ), (20) 
7*5 = K-<)/» + 0(Ax 2 ), (21) 

where e™, f£ are given by their continuum expressions fllOflll) . and the correction term is 
proportional to the derivatives of the fields and to the squares of the lengths of the cell 
edges. If the fields are slowly varying, meaning that the derivatives are small, the correction 
term can be neglected. This is what we mean by the continuum limit. 

We also need the discretized version of the curvature tensor T^: it is a plaquette. In 
our case the plaquettes are triangles, and we define the parallel transporter along a closed 
triangle spanning the k vertices: 

Pijk = UijUjkUki, Pjfjf = -t-Tt {jl AB Pijk) ■ (22) 

Expanding P^k around the center of the cell x we obtain: 

P ijk = 1 - - x»){xl - + 0(Ax 3 ) 

and 

P$S = ~ OK - x»)^ B (x) + 0(Ax 3 ). 

Using the above ingredients one can easily construct the lattice regularized version of the 
action terms (|T3l) . For example, the discretized cosmological term So has the form: 

e «ou---«d gAiA 2 ..A d 
all cells ^ ' ' 

where any number of e's can be replaced by /'s. In the continuum limit one uses Eqs.( l20ll2Tl) 
and obtains 

all cells ^ >' ' ' 

The coordinate factors combine into the volume of the cell (fT6l) and one gets 

S = V{cell) det(e) [1 + 0{Ax)) / rf d a; det(e) = 5 , (25) 

cells 

where det(e) is composed from the continuum tetrad ( JTUl) and is attributed to the center 
of a cell. Eq. (|25|) proves that the lattice action (|23|) becomes the needed continuum action 
f fT3|) if the fields involved are slowly varying from one lattice vertex to the neighbor ones. 
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Similarly, one finds the lattice version of all other action terms Sk of Eq. ([TBI) : 

EJ i,..i d f:A\A2...A d , v , N 

_ / pAiA 2 pA 3 A 4 pA 2 fc-iA 2fc \ / ~A 2k+1 ~A d \ n 

(d+l)\ d\ V * 0,ili2 io * 3i4 iohk-ihk J \ iahk+l ' ' hid) fc ' 

cells ^ 

(26) 

where the total number of plaquette factors P f[2"2"j) is fc, = 0, 1 . . . [d/2]. In fact one can 
write a variety of such action terms replacing any number of composite frame fields e f ll8p 
by the composite frame fields / (Tl9|) . 



C. Lattice partition function 

The lattice-regularized partition function for the spinor quantum gravity is quite similar 
to that of the common lattice gauge theory. One integrates with the Haar measure over link 
variables Uij living on lattice edges, and over anticommuting fermion variables living 
on lattice sites. The lattice, though, must be simplicial, otherwise the trick used e.g. in 
Eq. f[!Hj) to get the diffeomorphism-invariant action in the continuum limit, would not work. 

Because of the requirement of diffeomorphism-invariance, the lattice action is quite differ- 
ent from those used in common lattice gauge theory. Typically one has many-fermion terms 
in the action. There are no action terms without fermions. One can write 3 action terms in 
2c? (all of them are 4- fermion), 6 terms in 3d (four are 6-fermion and two are 2-fermion), 8 
terms in 4d (five are 8-fermion and three are 4-fermion), etc. We assume that spinor fields 
are dimensionless since we normalize the basic Berezin integrals as 

di/;^ = l, [d^^ = l, dip = 0, d^ = 0, (27) 



hence all quantities in Eq. ( 1261) are dimensionless. Therefore, the "coupling constants" 
one puts as arbitrary coefficients in front of the action terms Sk (I2"B"1) are all dimensionless. 
The partition function is 

Z= II JwUh II JdU ij exp^\Pst\^,^U)\ (28) 

vertices i links ij \ cells / 

where are lattice actions of the type (|26|) with any number of composite frame fields e 
( TIB]) replaced by the other composite frame fields / ( fT9|) . 
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V. MEAN-FIELD APPROXIMATION 



The partition function (128]) defines a new type of a theory, and new methods - exact, 
numerical and approximate - have to be developed. 

In principle, in order to compute the partition function (128 p as well as correlation func- 
tions, etc., one has to Taylor-expand the exponent in Eq. ( |28|) to certain powers of the 
fermionic action terms such that at all lattice sites there is precisely the same number of 
fermion operators ip^ and ip as there are integrations, since all other contributions are identi- 
cally zero by the Berezin integration rule (T27|) for anticommuting variables. The subsequent 
integration over link variables with the Haar measure is simple [7J since link matrices Uij 
never appear in a large power. Moreover, the majority of potentially possible contributions 
are killed by link integration. 

In practice, however, the arising combinatorial problem is tremendous, and we did not 
manage yet to find a computational algorithm that would be faster than the exponent of 
the lattice volume. So far we have done a toy model in Id exactly (see the Appendix) and 
succeeded in computing numerically correlation functions in 2d for limited volumes. There 
is a hope that the 2d model may be solved exactly but the method can hardly be extended 
to higher dimensions. 

Therefore, for this pilot study, we have developed an approximate mean-field method to 
get the first glance on the dynamics of the new interesting theory at hand. Comparing the 
results with an exactly solvable model we see that the mean field accuracy is within a few 
percent. In the 2d model there are a few exact functional relations that are satisfied with 
the accuracy better than 15%, and this can be systematically improved. More important, 
the mean field approximation reveals a nontrivial phase structure of the theory in the space 
of the coupling constants A^. This is the main finding of this study that may have important 
physical implications, see the Introduction. 

The mean-field approximation we use is an extension of methods developed in condensed 
matter physics, that go under the name of "dynamical mean-field approach" or "local impu- 
rity self-consistent approximation" or "cavity method" , see Ref. 2l| for a review. Roughly 
speaking, the idea of the method is the following: One first picks up a simple element of the 
lattice (e.g. one simplex or a group of simplices with or without the boundary, let us call 
this fixed element "the cavity"), and calculates the effective action for the fields inside the 
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cavity in the collective background of the fields outside it, replacing the background by the 
supposed mean field. At the second stage, one makes the method self-consistent, namely 
one calculates the mean field by integrating over the "live" variables inside the chosen cavity 
using the effective action found and expressed through the mean field at the first stage. As 
a result one gets a system of highly nonlinear self-consistent equations for a set of mean 
values of the field operators. Solving those equations one obtains the mean-field values as 
function of the coupling constants A^ . This gives the phase diagram of the theory in the 
space of the coupling constants. 

The method has the advantage that it can be systematically improved by enlarging the 
chosen cavity. In the limit when the cavity covers the whole lattice it is an exact calculation. 
Also, it is known that the accuracy of the mean field method is better the more nearest 



neighbors there are |2l|. In simplicial lattices the number of neighbor cells is large, and the 
mean field method becomes exact in the limit d — » oo. 

Let us formulate the method more mathematically. We choose the cavity, for example the 
elementary simplex with the boundary. We label the "live" fields belonging to the cavity by 
m,n, . . . and the fields outside the cavity (that will be replaced by mean fields) by . . .. 
The full partition function can be written symbolically as 

Z = I d^ m d^ m dU mn e Smn I dU mi e s ™ I d^jd^idUij e s ^ (29) 



where S mn is the part of the action that contains only fields from the cavity, Sij contains only 
fields from outside the cavity and S m i contains both. The link elements U m i are connecting 
vertices from the cavity with their nearest neighbors outside. 

The last integral in Eq. (|29l) is the full partition function with the cavity cut out. When 
the lattice volume goes to infinity cutting out a finite cell does not change the averages of 
operators as compared to the averages computed on a full lattice; we denote them as 

f dipjdipidUij e s » 

The integration over the links U m i connecting the cavity with the outside neighborhood 
must be performed explicitly in Eq. (T2"9"|) . We expand e Smi in powers of the mixed action; 
since S m i is a fermion operator the power series is finite. Integrating over U m i splits all terms 
involved into a sum of products of operators composed of the cavity fields O (ipln, if) m , U mn ) 
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and those living outside the cavity O' yfij, ipi, Ui 3 

I dU mi e s ™ = 1 + °p ^m, U mn ) 0' p tyl (31) 
J p 

where the sum goes over various fermion operators labeled by p. Operators built from 
the cavity fields are left intact whereas the outside operators are replaced by the averages 
according to Eq. ( )30l) . We, thus, obtain the effective action for the fields inside the cavity: 

e 5 cff , mn = e s mn (^1 + J2°p K» Vv», U mn ) ((OP) J . (32) 

Finally, we make the calculation self-consistent by requesting that the operator averages 
(O p ) computed from the cavity fields alone with the effective action (132]) coincide with the 
full ones ((O p )): 

ln v _ / djjldjj m dU mn e s ^^O p ip m , U mn ) 

^ P ' ~ r A /t A I ATT 9^ ~~ \\ U P/)- \ 66 ) 

J dY m d%p m dU mn e b ^,™ 

Since S e s depends on the averages ((O p )) the self-consistency Eq. (!33l is in fact a set 
of nonlinear equations on the mean values of the operators introduced in this derivation. 
Solving those equations one finds the values of the average operators as function of the 
coupling constants of the theory. 

Of special interest are the cases where certain operator averages (the "condensates") 
violate the continuous symmetries of the original theory. It signals the spontaneous breaking 
of symmetry and leads to a nontrivial phase diagram for the theory. In the next section we 
illustrate it in a general 2d model. 

VI. TWO-DIMENSIONAL SPINOR GRAVITY 

The partition function is defined by Eq. ( ]28|) where the action has in general three terms 
with three arbitrary coupling constants A 1,2,3, 



S = J rf 2 x^A 1 det(e) + A 2 det(/) + A 3 ^e^e^e^, A,B = 1,2. 



(34) 



The lattice- regularized version of it is, according to Eq. ( 123]) . 

Jjk P AB 



cells 
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where k = 0, 1, 2 label the vertices of a cell which in 2c? is a triangle. Using integration 
by parts the first term in (|34j) can be rewritten as — J^f (-^VO 2 - ft gives an alternative 
discretization for the same continuum action: 

ijk AJ3 / ' \ 

^ = EVV (-| Al ^ f( ^ <)2+Aa/ « / * +As ^ / *)- (36) 

cells ^ 

Although in the continuum limit the lattice actions (135]) and fl36|) differ by a full derivative, 
the lattice mean-field approximation gives numerically slightly different results depending 
on whether we start from Eq. (1551) or from Eq. (13"tTj) . The deviation serves as one of the 
checks of the accuracy of the approximation, and we find it consistent with other accuracy 
checks. 

In d = 2 the Lorentz group is the Abelian 5*0(2) ~ U(l) group. The spinors are two- 
component, and the 7- matrices are the Pauli matrices j A = a A , A = 1,2. The Lorentz 
rotations generator is S 12 = — o~ 3 /2, see Eq. (jSj). The analog of the "gamma-five" matrix in 
2d is 7 3 = ij lr y 2 = —a 3 . 

Both variants of the frame field e A and f A as well as there lattice extensions, and f^, 
are invariant under two global U(l)y x U(1)a transformations: 

vector transformation : ip — > e^tp, ip^ — > i/j^e~ 1 ^ , (37) 

■ cz 3 x x • a 3 

axial transformation : ip — > e l ^ a ip, ijy — > vjj'e 1 ?" . (38) 

Therefore, both the continuum ( 1341) and the lattice ( l35l) actions possess these two global 
symmetries also; the corresponding Nother currents are conserved. The vector symmetry 
means that the fermion number is conserved whereas the axial means that the difference 
between the numbers of "left-handed" and "right-handed" fermions (described by the upper 
and lower components of the spinors, respectively) is also conserved. It is also called the 
helicity conservation, or chiral symmetry. 

A. Exact results 

In the 2d partition function (1251) there are four integrals per site over fermion variables 
ip 1 , ip 2 , -025 an d one integration per link over the Abelian matrix U{j = exp (— i^-a 3 ) . 
Berezin's integrals over fermions (I27p are non-zero only when every lattice site takes exactly 
four fermion operators from the action exponent. Meanwhile, each term in the action ( l35|) 



17 



or (136]) is four-fermion. From counting the number of fermion fields coming from the action 
(which must be equal to the number of integrations) we conclude that the partition function 
Z is a homogenous polynomial of the coupling constants Ai.2,3 of order N, where N is the 
total number of sites in the lattice, 



C\\ x , A2,3 0, 

C 2 \?, Ai >3 -> 0, (39) 

, Al,2 — > 0. 



Since there are two types of frame fields, e and /, we can define three types of "physical" or 
invariant volumes of the generally curved space, averaged over quantum fluctuations of the 
composite frame fields, 

d , f ,2 < ■ ' - I dZ dlogZ 



(V,) = ( / d 2 xdet(e)) = ^ Wi = (40) 



d , I , 2 ■ I dZ dlogZ 



(V 2 ) = ( / d'xdet(f)) = ^Qx 2 = ^A^' (41) 



(Vs) = 



The immediate conclusion from Eqs.( l39ff4"2]) is that the average action is 

(S) = Xi(Vi) + X 2 (V 2 ) + \ 3 (V 3 ) = N = ^- (43) 

irrespectively of the coupling constants, where M is the number of triangle cells, which is 
twice the number of vertices N for large simplicial lattices. 

Further on, one can introduce "physical volume susceptibility" or variance 



(AV?) = {{V, - (V x )f) = (VI) - (V.Y = (44) 

(Avi) i «y 2 - iy 2 ) f) = (v?) - (v 2 f = *||* (45) 

(af 3 2 ) i ((y 3 - <^» 2 > = (v 3 2 ) - <^> 2 = (46) 

Therefore from Eq. ( 1391) we know exactly the average physical volumes and volume suscep- 
tibilities at least at the edges of the parameter space A = (Ai, A 2 , A 3 ): 

M M M 

\ v V\2,3^o = (^)a 1>3 -+o = (^3)ai, 2 ^o = (47) 



(AV?) M ^ = --^ (AV 2 2 ) Al ^ = (AV 3 2 ) Al ^ = (48) 
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where M = 2N is the total number of simplicial cells in the lattice. The proportionality of 
these quantities to M is a very general property (valid not only at the edges of the parameter 
space but everywhere) following from Eq. ( I39p . It shows that the physical volume is an ex- 
tensive quantity, as it should be. This is not altogether trivial since nonperturbative metric 
fluctuations allow, in principle, "crumpling" of the space, or the formation of "branched 
polymers", and that is what some researchers observe in alternative nonperturbative ap- 
proaches to gravity. In spinor gravity, it is a natural result following physically from the 
non-compressibility of fermions and mathematically expressed by Eq. (|39|) . 

The susceptibilities (jiSj) are also extensive, as should be expected. In the classical ground 
state there are no quantum fluctuations, so AV = 0. The fact that (1481) is nonzero means 
that we are dealing with a fluctuating quantum vacuum. At the same time for large volumes 
the relative strength of the fluctuations die out: V AV 2 /V ~ 1/a/M — > 0. 

There are theorems for mixed derivatives, valid in the whole parameter space A, that can 
be used to check the accuracy of approximate calculations, for example, 



Finally, there is an exact statement about the average curvature. The number of link 
variables in all terms of the action ( 135]) is even. That gives a nonzero result from integration 
over links for the partition function. However, if one attempts to compute the average of the 
curvature proportional to J 7 that in the lattice formulation is given by a product of three 
links (see Eq. (I22p ). the number of link variables becomes odd, and link integration yields 
an identical zero. Therefore, we conclude that the average Cartan curvature proportional 
to the average scalar curvature is zero, 



This result in 2d is, of course, in conformity with the zero Euler characteristic of a torus; 
no other result could be correct. It is illuminating, however, to see how "microscopically" 
the Euler theorem works for fluctuating spaces. In higher dimensions det(e)i? is not a full 
derivative but it still may be possible to find its average in a similar way. 

At the same time, owing to quantum fluctuations the average curvature squared is gen- 
erally nonzero and extensive, 



d (J d 2 xdet{e)) _ dlogZ _ d (J d 2 xdet{f)) 
d\ 2 <9Ai<9A 2 d\i 



(49) 



(det(e)E) = 2(^1) = 0. 



(50) 




M, 



(51) 
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implying that the volume- independent combination dies out in the thermodynamic limit, 



{(fd 2 xdet(e)R)) 1 

->0. (52) 



(/ d 2 x det(e)) 

Eqs.f l47|l50|52|) mean that although we apparently deal with a quantum fluctuating vacuum, 
the space is on the average large and flat in the absence of external sources. Therefore, one 
can say that the model describes a flat background metric G^ v that is a unity matrix in a 
particular frame representing the flat space but transforms as a tensor under the change of 
coordinates. We shall use this notion in Section VII. 



B. Mean-field approximation for one simplex cavity 

In this Subsection we apply the mean-field method formulated in Section V to the lattice 
action f )36|) where we first put for simplicity A 3 = 0. At the end of this Section we formulate 
the main results for A3 7^ 0. 

In the first approximation to the mean-field method, we choose the elementary triangle 
cell (m,n,p) as the "cavity", see Fig. 3. The fields inside the triangle cavity are consid- 
ered as real quantum fields, whereas the fields outside the cavity are combined into certain 
gauge-invariant operators that are frozen to their mean-field values. The triangle cavity is 
surrounded by three "black" triangles of the type (i,m,n) with a common edge, and by nine 
"white" triangles of the type (i,j,m) with a common vertex. The effective action for the 
fields inside the cavity gets contributions from both types of neighbors. 




FIG. 3: The simplest triangle cavity (m,n,p) and its neighbors used in the mean field calculation. 

Following the method of Section V, we expand the action exponent for every border cell, 
and integrate over the link variables U m { connecting the cavity with the outer lattice. As the 
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result we obtain the product of operators built of the fields inside the cavity and those built 
of the outside fields. The latter operators are replaced by the averages to be found later 
from the self-consistency condition. We stress that after integrating over Ui m the operators 
on both sides can be only gauge- invariant. 

For the single cell cavity we obtain operators of two types: single-site operators (they arise 
from the "white" cells), and double-site operators built from fermions at adjacent vertices 
(they arise from "black" cells). 

Here is the list of operators that appear in this calculation. First of all, there are operators 
that are invariant under the U(l)y x U(1)a transformations described in Subsection A: 



Old) 

o 2 (i,j) 



3/i 



(53) 



0\ is a single-site operator while O2 and O3 are double-site operators. 

To be able to study the potential breaking of the U(l)y x U(1)a symmetries we introduce 
operators that transform under those rotations. The chiral non-invariant operators that 
transform under U(1)a (jSHJ) are 

= i(4^), c l (i) = i{4a^ i ,) 

C 2 (i,j) = i [(Vj^) 2 (^Vi) + (^i) 2 (^M) 

c 2 (ij) = i [(^i) a (^Vi) + (^j) 2 (4^i) 

Fermion number violating operators transforming under U(l)y ( )37|) are 

Wi(i) = ipi,\ipi,2, 



(54) 



(55) 



All operators are Hermitian. 

The effective action for the fields inside the triangle cavity is computed as described in 
Section V. From the "black" neighbors we obtain the double-site effective action 



S{m,n) 
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A 2 [(0 1 (m) + 1 H)(0 1 ) + 3 ] 



+ 



Af 
36 



02(d) - 2 \Wl{Wi) + W 2 (Wt) 
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C2{C l ) + C 2 (C 1 ) 



(56) 



where all double-site operators refer to the cavity vertices m and n. From the "white" 
neighbors we obtain the one-site effective action 

e s{m) = 1 + |a?[20 1 (0 1 > + (0 3 >] 



\ 2 2 

H 

36 



1 (0 2 ) - 2 (W}(W 2 ) + W^Wl)) - Cx(C 2 ) + Ci(C 2 )] • (57) 



Actually the operator averages in ( 15B|57|) imply averaging over the cyclic permutation of 
lattice sites in the cavity: e.g. (Ox) = |(Oi(m) + Ox(n) + 0\(p)), and similarly for the 
double-site operators. The full effective action for the cavity is a sum over all 12 neighbor 
cells, 

gSeft = exp [s( m? n? p ) + n ) + S ( U} p ) + S '(p ) m )) + 3 (s( m ) + + 5(p))] ) (58) 



where S{m,n } p) is the original action for the cavity triangle (m,n,p), as given by Eq. 

We see that the effective action for the fields living in the cavity cell depend explicitly on 
the yet unknown operator averages (O), (C), (W). To find them, one equates the operator 
averages as defined by the effective action fl58l) to those introduced previously, see Eq. fl33l) . 
As the result one obtains a system of nonlinear self-consistency equations on the averages 
(O), (C), (W). Solving those equations one finds the averages as function of the coupling 
constants Ai >2 . 

This calculation is straightforward but the equations are rather lengthy. Therefore, we 
just comment here on its most important features. 

First of all, we notice that S e s is quadratic in the symmetry breaking operators Cx, 2 and 
Wx t 2, as it should be, therefore one gets a system of linear homogeneous self-consistency 
equations on the averages (Cx, 2 ), (Wx, 2 ) that always have a zero solution, unless the deter- 
minant of this set of linear equations is zero. If the determinant is nonzero in the whole 
range of the parameter space A there is no spontaneous symmetry breaking. If the determi- 
nant passes through zero at some surface in the A space, it is where the second order phase 
transition takes place. Inside the domain where one of the U(l) symmetries is spontaneously 
broken the condensates (Ci j2 ) or (Wi^) are nonzero and are found as anomalous solutions 
of the nonlinear equations, together with the symmetry-preserving averages (01,2,3). 

In the absence of symmetry breaking one puts (Cx t2 ) = (Wi, 2 ) = and solves the system 
of three nonlinear equations on the averages (01,2,3). There are in general several solutions 
but none are real in the whole A space. We pick up the solution that is real near the line 
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A2 = 0. However it develops a cut and becomes complex at the lines | A2 1 = 8.69 1 Ai | signalling 
that there can be a phase transition along these lines. A careful study in the next Subsection 
shows that indeed these are the border lines separating the phase with spontaneous chiral 
symmetry breaking, see Fig. 4. 
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FIG. 4: The phase diagram of the 2d spinor gravity in the (Ai,A2) plane at A3 = 0. Region / 
corresponds to the chiral-symmetry broken phase, region II is a regular phase. The dots show the 
lines of the 2 nd order phase transition: IA2I — 8.69|Ai|. 

Outside this domain, i.e. at | A2 1 < 8.69 |Ai| chiral symmetry is not broken, the solution 
for the normal, symmetry-preserving operators is real, and one can approach the line A2 = 
where we can check the accuracy of the mean field method by comparing the average physical 
volume (Vi) ~ X]ceiis(^ e ^(^)) wnere the average is computed over one (cavity) cell with the 
effective action fl58|) . with the exact result fj4"7"|) . We find numerically 

M M 
(^i)a 2 3->o — 0.572 — (mean field) vs. 0.5— (exact), (59) 
Ai Ai 

where M is the total number of lattice cells. We note that the functional dependence on 
Ai is correct whereas the numerical coefficient deviates from the exact one by 15%. A 
more powerful check comes from computing the average action (S) which turns out to be 
a constant up to the third digit in the whole range of analyticity in Ai 5 2,3, equal to 0.57, 
instead of the exact result ( 143]) being 0.5. This is the typical accuracy with which other 
checks with exact results are fulfilled. 

We have also tested a more primitive mean-field approximation where the cavity is taken 
in the form of two neighbor vertices connected by a link. It is also capable of detecting the 
spontaneous breaking of chiral symmetry but the accuracy is, of course, worse: it is at level 
of 40%. 
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C. Spontaneous chiral symmetry breaking 

An accurate way to study spontaneous breaking of a continuous symmetry is to introduce 
a small term in the action that violates the symmetry in question explicitly. Since we are 
interested in the spontaneous breaking of the U(1)a or chiral symmetry we introduce the 
simplest diffeomorphism-invariant "mass term" 



that is not invariant under chiral rotations (138|) . Its discretized lattice version is obvious, 
see the first term in Eq. ([35]) • 

Adding this term we repeat the same mean-field derivation of the effective action for the 
triangle cavity as in Eqs.f l56|l57|) which now obtain an addition 




Let us note that the terms linear in the mass parameter to are also linear in the chirality-odd 
operators C\^- 

With this addition to the previous effective action ( 1581) we now turn to solving the self- 
consistency equations for the operator averages (0^2,3) and (Ci )2 ), (^1,2) • At m 7^ there 
is a solution for the chiral condensates (Ci^) in the whole (Ai,A2) plane (we still keep 
for simplicity A3 = 0). However, the dependence of the chiral condensates on the mass 
parameter to is totally different depending on whether we are in the region I where chiral 
symmetry is broken, or in the region II where it is preserved. 

In the region II the dependence of the chiral condensates on the mass is linear at small 
to; if to goes to zero the chiral condensates vanish. In the region I the dependence of 
the chiral condensates on the small parameter to that breaks the symmetry explicitly is 
non- analytic. Actually the chiral condensates are proportional to the sign functions of to, 
(£1,2) ~ sign (to), see Fig. 5. The behavior of the chiral condensate (C1.2) in the whole 
(Ai,A 2 ) plane (at A 3 = 0) is shown in Fig. 6. 

Figs. 5,6 clearly demonstrate that there is a range of coupling constants where the theory 
undergoes spontaneous breaking of the continuous U(1)a chiral symmetry, with a line of the 
2 nd order phase transition separating the phases. 




(60) 





(61) 



(62) 
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FIG. 5: The dependence of the chiral con- 
densate (C\) on the mass parameter m with 
the varying value of Ai at fixed value of A2- 
The step-like behavior ~ sign(m) signals the 
spontaneous breaking of symmetry. In this 
example, Ai = 0.12 is the 2 nd order phase 
transition point where the chiral condensate 
abruptly vanishes when m = 0. 




FIG. 6: The value of the chiral condensate 
(C\) in the (Ai, A2) plane. At the phase tran- 
sition line I A2 1 = 8.69 I Ai| the condensate 
vanishes with an infinite first derivative. 



D. No fermion number violation 



The effective action ( 1581) contain operators W\$ violating the U(l)v symmetry (1371) re- 
lated to the fermion number conservation. If (Wi 2) 7^ it signals the spontaneous violation 
of this symmetry. Fermion number conservation is spontaneously broken e.g. in ordinary 
superconductors and "color" superconductors in QCD. However, in contrast to chiral sym- 
metry that is also broken in QCD, spontaneous fermion condensation usually happens not in 
the vacuum but at nonzero chemical potential for fermions, since interactions are effectively 
amplified near the Fermi surface. In this Subsection we look for the spontaneous fermion 
number non-conservation in the same 2d model where we observe spontaneous breaking of 
the U(1)a symmetry in the mean-field approximation. 

Following the same logic as in the previous Subsection we introduce a small action term 
that violates the U(l)v symmetry explicitly, 

SB-odd = J d 2 x det(e) bip 1 ip 2 . 

This operator preserves the chiral U(1)a symmetry. The correction to the effective one- 
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triangle action is: 



SB-odd n) 



SB-odd(m) 



A26 
27 

27 



(W 2 (0 1 )+0 3 (W 1 )), 



(0 1 (W 2 ) + W 1 (0 3 )). 



We solve again the self-consistency equations on the operator averages but now with this 
addition, and look for non-analytic dependence on the small parameter b. In contrast to the 
case of spontaneous chiral symmetry breaking we do not find such solutions in the whole 
parameter space A = (Ai, A2, A3). 

We conclude that the fermion number conservation is not broken spontaneously in the 
model, except maybe along the line of the chiral phase transition. We did not study the 
inclusion of a chemical potential for fundamental fermions - that would explicitly violate 
Lorentz symmetry but presumably make the phase diagram of the model more rich. 

There are no reasons why fermion number conservation would not break spontaneously, 
say, in Ad, and the mean field method suggested here is a simple way to detect it. 

E. Full phase diagram 

The full action compatible with the principles proclaimed has in 2d three terms and 
consequently three coupling constants. In the previous Subsections we have restricted our 
study to the case of A3 = 0. 

Actually we repeat all the steps described above also for A3 7^ 0. The algebra becomes 
more cumbersome but still doable. We find that the chiral symmetry breaking phase I 
occupies the cone 



shown in Fig. 7; Fig. 4 is its section at A3 = 0. We remark that the accuracy of the mean- 
field approximation for some reason deteriorates as A3 grows. Still the exact relation ( I43jl 
holds even at A 3 — > 00 up to a factor of 1.6. 

VII. LOW-ENERGY ACTION FOR PROPAGATING FIELDS 

The theory defined by the partition function (12*81) is in fact ultra-local: all correlation 
functions of gauge- invariant operators generally decay exponentially at the separation of a 



\\ < 77.23 \\ + 5.36 A3 



(63) 
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FIG. 7: Spontaneous chiral symmetry breaking takes place inside a cone in the full parameter 
space A = (Ai, A 2 , A 3 ). 

few lattice cells. This is clear on general grounds but we have also checked it by numerical 
simulations on a 2d lattice of limited volumes. Special measures should be taken to ensure 
that certain degrees of freedom propagate to distances that are large in lattice units. The 
situation here is different from the common lattice gauge theory where it is sufficient to take 
the limit [3 — > oo where [3 is the inverse gauge coupling to guarantee long-range correlations 
in lattice units. In our theory there is no such obvious handle. 

However there are ways to guarantee that long-range correlations appear; moreover that 
can be checked in the mean-field approximation. An example which we consider here is pro- 
vided by the Goldstone theorem: If a global continuous symmetry is broken spontaneously 
the associated Goldstone bosons are exactly massless and hence propagate to large distances. 

In the previous Section we have shown that the continuous U(1)a or chiral symmetry is 
spontaneously broken in a broad range of the space of the coupling constants. Supposing 
the coupling constants are chosen inside that range (inside the cone in Fig. 7), there is a 
massless Goldstone excitation a(x) being the phase of the U(1)a rotation ( l38i) . 

Under this rotation the chirality-violating operators Cx,i and transform as 

Cf = C x ± C t ->• e ±ia C±, Cf = C 2 ± C 2 -> e ±ia Cf. (64) 

To derive the low-energy action for the Goldstone field we allow the phase a to vary slowly 
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from cell to cell: 

(Ci ± 2 )=Pi,2e ±lQ(cell) . (65) 

We parameterize the operator averages (Cf 2 ) in the same way and re-derive the effective 
action (|58|) for the fields inside the triangle cavity, taking now into account that the operator 
averages have slightly different phases in the cells surrounding the cavity. Then, integrating 
over the fields inside the cavity we find the effective one-cell partition function Z\ modified 
by the varying nearest neighborhood. If a is the same for all neighboring cells it is the same 
expression as in Section V, let us call it Z\q. However, there will be further terms depending 
on the gradients of a(x), we are now after. 

The full partition function is, in the mean field approximation, a product of Z\§ over all 
cells whose number is M. Therefore the action for the Goldstone field a(x) is 



S G = -M\nZ l = -M In Zio (66) 
- M 



1 dZi If d 2 Z 1 I dZ x dZ\ \ , n , 3n 



Z w dai 2Z W \daidctj Z w dai daj 

where «j is the value of the phase attributed to one of the 12 neighbor cells i, and Aa, is the 
difference between «j and «o attributed to the central cavity cell; the summation goes over 
all neighbor cells. It is important that the dependence of Zi on ctj starts from quadratic 
terms, which is the consequence of chiral symmetry; hence dZi/dat = 0, and we are left 
with second derivatives. 

Ignoring the first a-independent term in Eq. fl66|) we find that the action is quadratic in 
the jumps Aa from one cell to the neighbor ones, 



S G = -M 1 ff Am Ac* j + 0(Aa 3 ). (67) 
2Z W daidaj 



We now introduce a coordinate system by mapping the centers of the cells to coordinates 
(Section IV). If the changes of a from a cell to neighbor cells are small we can expand 

Aa { = d^aAx 1 * + ^ d^aAx^ Ax\ + ... (68) 

where Axf = x^{i) — x M (0) is the distance between the coordinate attributed to the cell i and 
that attributed to the cavity cell, in a given coordinate frame x M (i). Putting this expansion 
into Eq. ( 1671) we obtain 

S G = -M ^J- AxfAx U j d,ad v a + 0(Ax 3 ). (69) 
2Z 10 daidaj J 
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The first factor M, the full number of cells on the lattice can be written as 

Ef d 2 x 
= I I 

cells 



M 



' V(ce\\) 

where V^cell) is the cell volume in a given frame, see Eq. ([1 

,. _J 1 d'Z, v 

a™o V(cell) Z 10 daidaj 1 j 



(70) 



The combination 



(71) 



transforms under the change of the map — > x'^{x) as a product of the contravariant tensor 
times the square root of the determinant of a covariant tensor, hence the notations in the 
right hand side or Eq. (J7TT) . Its particular form depends, of course, on the coordinate system 
chosen. For a concrete map to the Cartesian coordinates of the lattice drawn in Fig. 3 we 
find that it is proportional to a unity tensor, 



GG 



[IV 



(72) 



regular lattice 

where the proportionality coefficient T(Ai j 2,3) is shown in Fig. 8; it is proportional to a 
combination of the moduli of the chiral condensates p\$, see Eq. (!65|) . This result is in 
conformity with the average flatness of the space found in Section VI. A. If one chooses 
another coordinate map s/GG^ changes accordingly. 




FIG. 8: The normalization factor T(Ai, A2, 0) in the low-energy effective chiral Lagrangian, Eq. (172j) . 



We thus arrive at a diffeomorphism-invariant low-energy action for the massless Goldstone 
field: 

(73) 



Sg = \J d 2 x VGG" U d^ad v a . 
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This field can propagate infinitely far in lattice units since its masslessness is guaranteed by 
the Goldstone theorem. 

To complete our study of the spontaneous chiral symmetry breaking we derive the analog 
of the Gell-Mann-Oakes-Renner relation for the pion mass in QCD, expressed through the 
quark masses. If chiral symmetry is broken explicitly by a small fermion mass term (Ib'Up the 
phase of the chiral condensate becomes a pseudo-Goldstone field with a mass proportional 
to the square root of the fermion mass m. 

Indeed, the addition of the mass term (1601) changes the effective one-cavity partition 
function: 

Zx -> Z l + mZ m + 0(m 2 ). (74) 

The Z m piece depends explicitly on the chiral condensate phase a introduced in Eq. (1651) . 
and from symmetry considerations it is clear that the expansion starts from the a 2 term; 
direct calculation confirms it. 

Summing up the mean-field action over the whole lattice one uses the relation (1701) where 
the r.h.s. behaves as ~ J d 2 x \/G according to the transformation properties under the 
change of the coordinate system x^(i) attributed to the lattice. We obtain thus the action 
for the pseudo-Goldstone mode in the continuum limit 

S G = X - J d 2 xVG (G^ d»ad v a + fx 2 a 2 ) + 0{d 2 ad 2 a) + £>(a 4 ), fi 2 ~ m, (75) 

where \i is proportional to the pseudo-Goldstone boson mass. We see that it is proportional 
to the square root of the mass parameter m that breaks chiral symmetry explicitly. In QCD, 
this is known as the Gell-Mann-Oakes-Renner relation for the pion mass. The coefficient 
in this relation depends on the coupling constants Ai 2,3- At the 2 nd order phase transition 
surface of the cone in Fig. 7 the pseudo-Goldstone mass goes to zero at fixed m. 

There is a famous Mermin- Wagner theorem stating that a continuous symmetry cannot 
be spontaneously broken in 2d as the resulting Goldstone bosons would have an unaccept- 
ably large, actually divergent free energy. Since the mean field approximation misses the 
Goldstone physics, one can argue that the spontaneous chiral symmetry breaking we observe 
is an artifact of the approximation. If, however, the Goldstone field a(x) is Abelian as here, 
the actual phase is, most likely, that of Berezinsky-Kosterlitz-Thouless where the chiral 
condensate pe ia indeed vanishes owing to the violent fluctuations of a(x) defined on a circle 
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(0, 27r), but the correlation functions of the type (e ia(x ) e~ ia ^) have a power-like behavior, 
and there is a phase transition depending on the original couplings of the theory. 

In any case our primary goal here is to learn how to deal with the lattice regularized 
spinor quantum gravity which is a new type of a theory. The mean-field approximation is 
one possible approach that is expected to work even better in higher dimensions where, as 
a matter of fact, the Mermin- Wagner theorem does not apply. 



VIII. HOW TO OBTAIN EINSTEIN'S LIMIT? 



The apparent diffeomorphism-invariance of Eq. (175]) is built in by our construction of the 
lattice and lattice action in Section IV. As soon as there are degrees of freedom that can 
propagate to long distances, their low-energy effective action is diffeomorphism-invariant in 
the continuum limit. 

In the previous Section the appearance of long-propagating mode has been guaranteed by 
the Goldstone theorem. However, it concerns only the specific Goldstone modes associated 
with the spontaneous breaking of continuous symmetry. Other degrees of freedom remain 
heavy: their correlation functions decay exponentially after a few lattice cells. If one at- 
tempts to write an effective low-energy action for the classical metric tensor g^ v (see below 
its exact definition) it will have the diffeomorphism-invariant form, 

Slow = J dx y/g*(- Cl + c 2 R(g cl ) + ...), (76) 

with the constants computable, in principle, from the original coupling constants of the 
lattice-regularized theory. However, if one does not take special measures, the ratio \ fc\j~C2, 
playing the role of the graviton mass, will be on the order of the inverse lattice spacing. In 
such a situation it is senseless to introduce the metrics in the first place. It makes sense only 



if \fc\fci happens to be zero or very small, such that the graviton and the Newton force 
propagates to large distances. 

To ensure it, it is sufficient to stay, e.g., at the phase transition of the second order, where 
all degrees of freedom become massless. The classical metric tensor g c ^ v and the effective 
action functional rfg^J can be introduced by means of the Legendre transform (proposed 
in this context also by Wetterich [14]). One introduces first the generating functional for 



31 



the stress-energy tensor Q^ v as an external source 



e 



j #t dt p ^ eX p (s + j g^ &A , (77) 



where S the fermionic action and g^ v is a 4-fermion operator built from the frame fields 
(I10|llll) or, after discretization, from their lattice versions (I18|I19I) . The classic metric field is 
by definition 

cl A SW[Q) 

This equation can be solved to give the functional 0^[g c1 ]. Using it one can construct the 
effective action as the Legendre transform: 

T[g d ]=W\Q}-g%^. (79) 

At the phase transition fluctuations are long-ranged. For long-range fluctuations it is legal to 
take the continuum limit of the lattice, which is diffeomorphism- invariant. The low-energy 
limit of diffeomorphism-invariant actions for a quantity transforming as a metric tensor is 
uniquely given by Eq. f J76|) . Moreover, the cosmological term necessarily has zero coefficient, 
Ci = 0, since otherwise the graviton would propagate to a finite distance \fcqjc\-, which 
contradicts the masslessness of the fluctuations at the phase transition. This is how one can 
recover Einstein's gravity from the lattice-regularized spinor theory. 

In principle, the effective Einstein-Hilbert action from spinor quantum gravity can be 
derived in the mean-field approximation similarly to our derivation of the low-energy effective 
chiral Lagrangian in Section VII. However, in 2d where we have so far succeeded in developing 
the mean-field method the Einstein-Hilbert action is a full derivative and there are no 
gravitons or the Newton force. Therefore, the derivation of Eq. (1761) has to be postponed 
till higher dimensions are studied along the lines of the present paper. 



IX. DIMENSIONS 



In this paper, we use unconventional dimensions of the fields, which, however, we be- 
lieve are natural and adequate for a microscopic theory of quantum gravity. The fermion 
fields are normalized by the Berezin integral (|271) and are dimensionless, hence the compos- 
ite frame field (fT0|) has the dimension 1 /length and the metric tensor has the dimension 
1/length 2 , in contrast to the conventional dimensionless metric tensor. On the other hand 
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all diffeomorphism-invariant quantities are dimensionless in our approach. In this section 
we explain why it is convenient, and what is the relation to the usual approach. 

The historic tradition in General Relativity is that the space-time at infinity is flat, there- 
fore one can safely choose the coordinate system such that is a unity matrix there. This 
sets the traditional dimensions of the fields. In particular, the scalar curvature has the 
dimension 1/length 2 , the fermion fields have the dimension 1/length^, etc. However, in a 
diffeomorphism-invariant quantum theory where one can perform arbitrary change of coordi- 
nates x M — » not necessarily identical at infinity, for example, a dilatation x^ — > x^/p, 
and where g^ v can a priori strongly fluctuate at infinity, this convention is not convenient. 

The natural dimensions of the fields are those that are in accordance with their trans- 
formation properties: any contravariant vector transforms as x^ 1 and has the dimension of 
length, a covariant vector, in particular the frame field transforms as a derivative and has 
the dimension 1/length, g^ u has the dimension 1/length 2 , etc. World scalars like the scalar 
curvature and the fermion fields are, naturally, dimensionless. In fact it is a tautology: a 
quantity invariant under diffeomorphisms is in particular invariant under dilatations and 
hence has to be dimensionless. 

In this convention, any diffeomorphism-invariant action term is by construction dimen- 
sionless and is accompanied by a dimensionless coupling constant, as in Eq. (|28|) . 

Let us suppose that we have a microscopic quantum gravity theory at hand that success- 
fully generates the first terms in the derivative expansion of the effective action, 

T = - Cl j d 4 x ^g + c 2 J d 4 x y/gR + ... (80) 

where Ci i2 are certain dimensionless constants expressed through the dimensionless couplings 
^i,2,... of the original microscopic theory. The ground state of the action (1) is the space 
with constant curvature R = 2ci/c2, represented e.g. by a conformal-flat metric 

^ = (((x-^ + p 2 )) ^ (81) 
where Xq and p are arbitrary. At the vicinity of some observation point Xq it can be made 
a unity matrix by rescaling the metric tensor, 



> _ r /6C2 2 

9v» = m g„ v , g„ v = <V, m = y — - (82) 

where the rescaling factor m has the dimension of mass, and g^ has the conventional zero 
dimension. At this point one can rescale other fields to conventional dimensions, in particular 
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introduce the new fermion field ip of conventional dimension m 3 / 2 : 

^ = m - 3/2 ^, ^=m- 3/2 ft. (83) 

The new composite dimensionless tetrad field compatible both with Eqs. ( 1821) and ( 1831) is 

= i 4 = ij i (V 7 A V^ + <FV M rV) . (84) 



One can now rewrite the action flHUj) together with the fermionic matter in terms of the new 
rescaled fields denoted by a bar, 

S = -cim^Jd A x^+ crm^ J d 4 x + m° J d 4 x y/§e A " VV M $ + h.c.) . 

2A=A 4 M|=l/V167rG N 

(85) 

Underbraced are the cosmological constant and the Plank mass squared, respectively; nu- 
merically, A = 2.39 ■ 1CT 3 eV, Mp = 1.72 ■ 10 18 GeV. The dimensionless ratio of these 
values, 

is the only meaningful quantity in pure gravity theory, independent of the arbitrary scale 
parameter m. If a fermion obtains an effective mass, e.g. as a result of the spontaneous 
chiral symmetry breaking, leading to an additional term in the effective low-energy action 

S m = J d 4 xy/g^ Mip = mM, J d 4 x VI W , (87) 

fermion mass m j 

then the "theory of everything" has to predict also other dimensionless ratios. For example, 
taking the top quark mass m t = 172 GeV one has to be able to explain the ratio 

m t M 



VAMp 



i i 



0.0848. (88) 



In other words, one can measure the Newton constant (or the Planck mass) or the cosmo- 
logical constant in units of the quark or lepton masses or the Bohr radius. Only dimensionless 
ratios make sense and can be, as a matter of principle, calculated from a microscopic theory. 
To that end it is convenient and legitimate to use natural dimensions when g^ u has the 
dimension 1/length 2 whereas all world scalars are dimensionless, be it the scalar curvature 
R, the interval ds, the fermion field ip or any diffeomorphism-invariant action term. 
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X. CONCLUSIONS 



We have formulated a lattice-regularized spinor quantum gravity that is well defined and 
well behaved both for large-amplitude and high-frequency fluctuations. In any number of 
dimensions one can construct a variety of fermionic actions that are invariant (i) under local 
Lorentz transformations and (ii) under diffeomorphisms in the continuum limit. We have 
built quite a few action terms satisfying (?) and (ii) for any number of dimensions. In fact 
our list of possible fermionic action terms can be expanded further if some of the additional 
requirements are relaxed. Therefore, we actually formulate a whole class of new kind of 
theories in any number of dimensions, characterized by a set of dimensionless coupling 
constants Ai i2 ,.... 

The continuum limit shows up if all degrees of freedom or at least some of them are slowly 
varying fields from one lattice cell to another. This is, generally, not fulfilled: generically, all 
correlation functions decay exponentially over a few lattice cells. For such "massive" degrees 
of freedom the theory is at the "strong coupling" regime where the continuum limit is not 
achieved and remains dormant. 

There must be special physical reasons for massless excitations in the theory, for which the 
continuum limit makes sense and diffeomorphism-invariance becomes manifest. One such 
reason is spontaneous breaking of continuous symmetry where the existence of massless 
fields is guaranteed by the Goldstone theorem. To show that spontaneous breaking may be 
typical in such kind of theories, we have developed a new mean-field approximation. We 
have checked its accuracy in a Id exactly solvable model, and in a full 2d theory where 
certain exact relations can be derived. The exact relations tell us nice things: the physical 
or invariant volume occupied by the system is extensive as due to the non-compressibility 
of fermions, the volume variance (or susceptibility) is also extensive showing that it is a 
true quantum vacuum, and the average curvature is, at least in 2d, zero meaning that the 
quantum space is on the average flat. They also tell us that our mean- field approximation 
is rather accurate, and the accuracy can be systematically improved. 

We show, within the mean-field method, that the spontaneous breaking of chiral symme- 
try happens in a broad range of the coupling constants and that in this range the low-energy 
action for the Goldstone field (or pseudo-Goldstone if we add a term explicitly breaking sym- 
metry) is diffeomorphism-invariant, as expected. 
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To obtain the low-energy Einstein limit one has to stay at the second-order phase tran- 
sition surface in the space of the coupling constants. There the masslessness of excitations, 
and not only of the Goldstone ones, is guaranteed. Hence one can go to the continuum limit 
where the diffeomorphism-invariance is also guaranteed by construction. Therefore, we ex- 
pect that the effective low-energy action for the classical metric tensor, derived through the 
Legendre transform, is just the Einstein-Hilbert action, with the zero cosmological term. 
This can be probably seen already in the mean-field approach for dimensions higher than 
two. This work is in progress. 
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Appendix. MEAN FIELD METHOD IN A Id MODEL 

We consider here a Id toy model with fermions and U(l) gauge symmetry, that can be 
solved exactly. We then apply the mean field method to this model to check how accurately 
does it reproduce the exact solution. We obtain quite satisfactory results. 

We take the same fields as in the full 2d model, namely the doublet of fermion fields 
ip a , iftlt, a = 1,2, that transform under the U(l) gauge transformation as ip ~ ► V* ~* 
ip^V\ V = exp(iaas). The gauge field is represented by link variables Z7y = exp(— iuiijcr 3 /4) 
that transform as Z7y — > ViUijVj . 

We construct the lattice version of the two "frame" fields, as in Eqs. (ll8|ll9j) . 

~4i = Wpn^Ui^j - #7^), A = 1,2, 

ft = 4^ A Ui^j - i^pji^i, (A.l) 

(which however do not have the meaning of frame fields in Id), and form the action that is 
quite similar to the full 2c? action (136 j) : 

N r 
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(A.2) 



where Ai^ and // are the coupling constants. The partition function is defined as a product 
of Berezin integrals on the Id lattice with N points: 

N 

Z = E[ / dipld^diplidipUdUij+i e s . (A.3) 
i=i 

We imply antiperiodic boundary conditions for fermion fields and periodic boundary condi- 
tions for link variables. 

The partition function (1A.3j) is exactly computable by a kind of transfer-matrix method. 
Diagonalizing the transfer matrix we obtain a nontrivial result: 

Z = 2(1 + (-1)*)A^ + (\ x - v / 2^/a? + A| + ^ + f Ai + v^A? + A| + ^ . (A.4) 

The fact that the partition function has the form of a sum of extensive exponents means that 
actually it describes simultaneously four independent phases or states of the system that do 
not compete and hence do not mix up in the thermodynamic limit A" — > oo. The stable 
phase is the one with the lowest free energy, that is with the largest partition function. 

Depending on the relation between the coupling constants one of the terms in Eq. (1A.4I) 
prevails at A" — > oo: 

• Phase 0: | A2 1 > |Ai| and — (Af + A 2 .) < fi < —\ [\\ ± A 2 ) 2 ; the partition function is 
given by the first term 

• Phase 1: /i > — | (Ai + A 2 ) 2 , Ai < 0; the partition function is given by the second term 

• Phase 2: /i > — | (Ai + A2) 2 , Ai > 0; the partition function is given by the third term 

• Phase 3: /i < — (A 2 + A 2 .), |Ai| > | A2 1 ; the partition function is complex and has no 
smooth thermodynamic limit. 

In phase fermions in the neighbor lattice sites form an ordered state of the type "pair-gap- 
pair-gap..." all over the lattice, where the "pair" means that there are four link matrices in 
the integration over link variables, and "gap" means zero matrices. It can be realized only 
on even- A" lattices, hence it is a lattice artifact, and we do not consider it further. Phases 1 
and 2 are states where two link matrices appear in all link integrations. We concentrate of 
phases 1 and 2 only in what follows. 
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We also calculate exactly average values of the following 4-fermion operators: 



4 — 7== = = m phase 1 

1 4 — X = = m phase 2 , 
V^I+M2a|+a|+2^) ^ ' 



. v / 2A 2 (A 1 + v / 2 v / A ?+ A 2+A t ) • r. i 

i 4 — / o o \ \ 111 phase 1 

(fAfA\ J VAfAl+/^(2A| +Aj+2 At ) /. fi x 

a "- ;/ 1 4 ^(-^v^I^) inphase2; 1 ' j 

A/A|A|+^A|+Af+27t) 



-1 



in phase 1 



— — = in phase 2 . 

^ 2 v / 2a/a|+a|+^ 

We now turn to constructing the mean field approximation to the model, to check its 
accuracy against the exact calculation. We apply the general method of Section V, which 
is rather straightforward in this simple case. We first take the "cavity" in the form of 
two neighbor lattice sites connected by a link (1 st approximation), and then three adjacent 
sites connected by two links (2 nd approximation). Both mean-field approximations give 
satisfactory accuracy when compared to the exact results but the second is, of course, better. 

In both cases the cavity boundary is just the neighbor sites connected to the cavity by link 
variables U m %- Expanding e Smi up to the second power (higher powers are zero because of 
too many fermion operators) and integrating over U m i we obtain several operator structures. 
Splitting them into the product of operators composed of fields inside the chosen cavity, and 
the operator built of the outside fields, we replace the latter by the averages, to be found 
self-consistently. Most of the operators break the <r 3 -symmetry of the original action, and 
we ignore them. The only operator with proper symmetries left is O = (tpl'ipi) 2 . We find the 
effective action for the two-site cavity 

e sM = 1 + Xl ( (0) + {^ m f) + 2 (A? + \\ + v)(0)(^ m ) 2 . (A.8) 

To obtain the self-consistency equation we equate the average of O found from the effective 
action (]A.8j) . to (O). The resulting nonlinear equation on (O) has three solutions. We 
choose the solutions (0)(Ai, A2, jj) that are real in the ranges 1 and 2 above. 

The results for the averages of three operators in the first and second approximations as 
well as their exact values (1A.5IIA.6IIA.7I1 are presented in Fig. 9. There is a "phase transition" 
between phases 1 and 2 at Ai = 0. We see that the second mean-field approximation 
corresponding to a three-site, two-segment cavity gives a very satisfactory accuracy. 
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FIG. 9: Averages of three operators in the Id model. The bold line is the exact result, the dashed 
line is the result of the 1 st mean- field approximation using one lattice segment, the solid line is the 
result of the 2 nd (two-segment) mean-field approximation. 
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